Surface-electronic structure of La(0001) and Lu(0001) 
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Most spectroscopic methods for studying the electronic structure of metal surfaces have the dis- 
advantage that either only occupied or only unoccupied states can be probed, and the signal is cut 
at the Fermi edge. This leads to significant uncertainties, when states are very close to the Fermi 
level. By performing low-temperature scanning tunneling spectroscopy and ab initio calculations, 
we study the surface-electronic structure of La(0001) and Lu(0001), and demonstrate that in this 
way detailed information on the surface-electronic structure very close to the Fermi energy can be 
derived with high accuracy. 

PACS numbers: 73.20.At, 71.20.Eh, 73.50.Gr, 68.37.Ef, 71.15.Mb 



After the first theoretical prediction of a Tamm-like 
surface state on Gd(0001) surface states with d z 2 
symmetry around the center of the surface Brillouin zone 
(BZ) have been found on all trivalent lanthanide met- 
als by photoemission (PE) and inverse photoemission 
(IPE) |1 H H H H 0. In all cases, the surface state 
was found to be located very close to the Fermi energy, 
Ep, where the accuracy of these experimental meth- 
ods deteriorates due to the signal cutoff at Ep. This 
weakness is manifested in the case of Gd(0001): First, 
a combined PE and IPE study concluded that the ex- 
change splitting of the surface state decreases with in- 
creasing temperature and vanishes within experimental 
accuracy above the Curie temperature (Tc) later, 
a scanning tunneling spectroscopy (STS) study showed 
that in fact a constant residual splitting persists above 
Tc 0. This was the first striking demonstration that 
STS is a very sensitive technique when states are located 
very close to Ep, because both occupied and unoccu- 
pied states can be probed within a single measurement, 
without any influence on the spectroscopic signal when 
crossing Ep. Consequently, recent systematic studies of 
magnetic lanthanide metals across the whole series of lan- 
thanide metals demonstrated that STS at low tempera- 
tures can probe these surface states with unprecedented 
accuracy pTol fllL . 

Here, we focus on the nonmagnetic 4/ metals La and 
Lu. In both cases, PE and IPE could not clarify whether 
the binding energy at the center of the BZ, T, is below or 
above Ep. While for Lu(0001) only PE data exists 0,111, 
a combined PE and IPE study of La(0001) concluded 
that the surface state is partially occupied 5j. However, 
due to the total system resolution in these experiments, 
the reliability of this latter statement must be consid- 
ered to be questionable. Furthermore, theoretical band 



structures are not available for La(0001) and Lu(0001) 
up to now. We therefore performed low-temperature STS 
studies and band-structure calculations for these two lan- 
thanide surfaces. In the present work, we compare the 
theoretical and experimental results, and we show that 

- with the help of the calculated surface-band structure 

- the STS data can be described accurately by simple 
models which do not only yield surface-state binding en- 
ergies but also lifetime widths and the dispersion of the 
surface bands. 

The experiments were performed in an ultrahigh vac- 
uum (UHV) chamber equipped with a low-temperature 
STM operated at 10 K [l(|- The samples were pre- 
pared in situ by vacuum deposition of electron-beam- 
evaporated La or Lu metal (« 30 ML), respectively, on 
a clean W(110) single crystal substrate followed by an- 
nealing at 800 K in case of La and 1000 K in case of 
Lu. STM images were taken in constant-current mode, 
and the STS spectra were recorded with fixed tip posi- 
tion, i.e. open feedback loop. The tunneling current / 
and the differential conductivity, dl/dV, were recorded 
as a function of sample bias voltage, V, by modulating 
V and recording the induced modulation of I via lock-in 
technique. A modulation amplitude of 1 mV (rms) at a 
frequency of w 360 Hz was used, with the time constant 
of the lock-in amplifier set between 10 and 100 ms, at a 
sweep rate of sa 2 mV/s. The spectra were taken in both 
directions, from lower to higher and from higher to lower 
sample bias, in order to verify that energy shifts (due to 
the finite time constant) are below 1 meV. Since both the 
STM tip and the sample were cooled to 10 K, the energy 
resolution was ~ 3 meV, corresponding to 3.5/cbT '. 

The calculations were performed using density func- 
tional theory (DFT) in the generalized gradient approx- 
imation (GGA) as given by Perdew and Wang [l4|. We 
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FIG. 1: (a) Tunneling spectrum of La(0001) at T = 10 K 
(rotated) and calculated band structure of a relaxed 11 layer 
La(0001) film. Surface states are marked by filled circles, 
bulk states by open circles. The most pronounced peak in 
the STS data at around 0.1 eV is attributed to the narrow 
surface-state band around the T point of the surface Brillouin 
zone, (b) Tunneling spectrum of Lu(0001) at T = 10 K and 
corresponding band structure of a relaxed 12 layer Lu(0001) 
film. The surface state appears as a sharp peak at Ef in STS. 



use the full-potential linearized augmented plane-wave 
method in film geom etry |l5l Il6| as implemented in 
the FLEUR code [IT]]- Spin-orbit coupling is induced 
self-consistently as described in Ref. ,18. For a proper 
description of the 4/ electrons, we apply the LDA+U 
method |l9l in an implementation similar to that of Shick 
et al. 20]. For La, we used values of U — 8.1 eV and 
J = 0.6 eV, for Lu of U = 4.8 eV and J = 0.95 eV. 
These parameters were chosen to simulate the experi- 
mentally observed positions of the 4/ bands |2l|. The 
La(0001) surface was simulated by an 11-layer film em- 
bedded in semi-infinite vacua, for the Lu(0001) surface 
a 12-layer film was used. The vacuum region was cho- 
sen to start 3.2 a.u. and 2.9 a.u., respectively, above the 
surface atoms of La and Lu. The muffin-tin spheres had 
a radius of 3.0 a.u. in case of La and 2.8 a.u. in case of 
Lu. A plane-wave cutoff of K max = 3.4 (a.u.) -1 was used, 
and the irreducible part of the two-dimensional BZ was 
sampled with 21 special fcii points. 

In Fig. [Ha), the STS spectrum of La(0001) is shown 
(rotated by 90°), together with the calculated band struc- 
ture. The "overview" spectrum is dominated by a narrow 
peak at V ~ 0.1V which is the signature of the <i z 2-like 



surface state. As STS probes mainly states around the 
r point of the projected surface BZ, the interpretation 
of the peak as a surface state is clearly confirmed by the 
band structure calculation, where a narrow downward- 
dispersing surface-state band is clearly visible in the local 
gap around T at about 0.2 eV above Ep. Fig. \Rb) dis- 
plays the STS data and corresponding calculated band 
structure of Lu(0001). Again, the only feature in STS is 
a sharp peak located close to V = (i.e. E — Ep)- It can 
be identified as the surface state by comparison with the 
band structure, which shows a surface state at 90 meV. 

A comparison of the band structures of La(0001) and 
Lu(0001) with those of other lanthanide metals shows 
that their electronic structures around T are very simi- 
lar, i.e. only small systematic changes occur across the 
lanthanide series. A clear trend in the band dispersion 
is visible: While for La, the effective mass is negative 
(i.e. the dispersion is downward towards lower energies 
E — Ep), Lu first shows a weak upward dispersion at 
r, which changes to downward with a band maximum 
at « 1/4 way towards the BZ boundary. The disper- 
sion of Gd is "in between" with a very flat surface-state 
band 123]. Calculations for Ce and Tm confirm this 
trend l23j. This systematic change further leads to 
a decrease of the surface-state bandwidth and higher ef- 
fective masses across the lanthanide series. From Fig. ^ 
\m*/m\ > 2 for La and \m*/m\ > 5 for Lu can be es- 
timated. Experimentally, in another STS study of Gd, 
Ho, and Lu we estimated \m* /m\ > 5 ^(j; a recent PE 
study of Ce concludes that \m*/m\ w 7.4 |2^]. These 
large effective masses are the result of the high degree of 
lateral localization of the Tamm-like lanthanide-surface 
states. Consequently, they appear as peaks in STS rather 
than as step-like functions (as would be expected for de- 
localized Shockley-like surface states, e.g. on the (111) 
surfaces of the noble metals). 

For a quantitative analysis of the spectral shape of the 
surface states in STS, we have applied a fit analysis that 
is based on planar tunneling theory |24|. The fit model 
was already introduced in Ref. ITol Therefore, we shall 
only briefly summarize it here. Assuming that the DOS 
of the tip is constant (within the energy range of inter- 
est) and that the transmission coefficient for tunneling, 
T(E), is bias-independent (for small V), the differential 
conductivity is given by 

Jj^cx / (n s T)(E)f'(E-eV)dE, (1) 

— oo 

where n s is the DOS of the sample, and /' is the 
derivative of the Fermi-distribution function; the latter 
takes into account the thermal broadening of the spec- 
tra due to the finite temperature in the experiment. 
For a surface band with quadratic dispersion, E{k\\) = 
Eo + (h 2 /2m*)fc|, and negative effective mass (m* < 0), 
the DOS is a step function: n s cx 9(Eq — E). Further- 
more, assuming |m*| m and E±_ = E — Eu < (j) Q g, T 
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FIG. 2: Highly resolved tunneling spectrum of La(0001) at 
T — 10 K and least-squares fit (grey line) composed of a 
surface-state line according to Eq. to (dashed) and a 
quadratic background (dash-dotted) to account for tunnel- 
ing into bulk states. The residual of the fit is shown as a thin 
dotted line. 



FIG. 3: Highly resolved tunneling spectrum of Lu(0001) at 
T — 10 K and least-squares fit (grey line) composed of a 
surface-state line according to Eq. to (dashed) and a 
quadratic background (dash-dotted) to account for tunnel- 
ing into bulk states. The residual of the fit is shown as a thin 
dotted line. 



can be approximated as: 

T{E) = exp(-2z^(0 cS ~E ± )) 

= cxp(-2z^(0 c ff-£) + fcj!) ( 2 ) 
oc exp [-pi(E - E)] , 

where cf> c s is the effective barrier height, z is the tip- 
sample distance, and E± is the perpendicular energy 
component. p\ is a constant parameter, with pi — 
zy / 2m/h 2 c ff (to* /to) [To| . Thus, the important term in 
Eq. ^ n s T, is simply an exponential function that is cut 
at the band maximum E by the step function. 

To account for finite lifetimes, each contributing state 
of the surface band is broadened by a Lorentzian func- 
tion. This is realized by writing n s T as an integral over 
Lorentzians weighted with T: 

n s T{E) oc / (B _ £) 2^ ( ( r (£)/2) 2 X ^ 
6(E - e) exp [-px(E Q - e)] de. 

r is the full width at half maximum of the Lorentzian, 
which is related to the lifetime r by V — H/t. T 
is written as a function of energy in order to account 
for possible energy dependences of the lifetime due to 
electron-electron (e-e ) and electron-phonon (e-ph) scat- 
tering [25l l26l I27L |2^ | . 

In summary, this model reflects the influence of band 
dispersion which leads to slightly asymmetric peaks with 
broad leading and narrow trailing edges. Although the 
above approximations are not valid for a strongly dis- 
persing band, this model can be qualitatively conveyed 
into the model of Li et al., who analyzed the width of 
the Ag(lll) surface-state spectrum [29j. While a large 
Pi factor (due to a high effective mass) leads to narrow 
peaks in STS, a small m* (as for the sp-like surface states 
of the noble metals) would result in a small p\. So the 



exponential decrease of the STS signal becomes negligi- 
ble, and a step function remains which is broadened by 
Lorentzians (corresponding to an arctan function, as in 
Ref. l2flh . Indeed, on a larger energy scale also in tun- 
neling spectra of the surface states on the (111) faces of 
noble metals a deviation from a step function can be ob- 
served in the form of a weak decrease of dl/dV towards 
higher energies [30ll3l"l|. 

In Fig. [3 a highly resolved tunneling spectrum is plot- 
ted together with the results of the fit analysis. As 
can be seen from the residual, the fit describes the ex- 
perimental curve very well. This confirms that - as 
for Gd, Tb, Dy, Ho, and Er pH d - the applied 
model works convincingly for a monotonically downward- 
dispersing surface state with high effective mass. Ac- 
cording to the fit result, the band maximum is located at 
Eq = (130 ± 5) meV, and the lifetime width at the band 
maximum is T(i?o) = (49 ± 10) meV, which corresponds 
to a lifetime of t(Eq) — (13 ± 3) fs. Thus, the band is 
clearly unoccupied at T and crosses the Fermi level at 
about 1/4 way towards the BZ boundary. The exper- 
imental value for the band maximum is about 70 meV 
below the theoretical one. As can be seen in the compar- 
ison with Ce(0001) [13, Gd(0001) H, and Lu(0001) (see 
below), DFT calculations of the unoccupied lanthanide- 
surface states always yield slightly higher energies for the 
band maximum. 

In the following discussion, we compare the results 
with literature. A combined PE and IPE study of 
La(0001) showed that the surface state is partially oc- 
cupied at T at room temperature 0, As this is 
clearly not the case in the present STS study at 10 K, 
the former observation can only be explained as a tem- 
perature effect in two ways: (i) The surface-state band 
may shift down with temperature. However, we have also 
applied STS on La(0001) at T w 60 K and did not see 
any shift of the band maximum. Therefore, we exclude 
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FIG. 4: (a) Same tunneling spectrum of Lu(0001) as in Fig.|3 
The grey line represents a least-squares fit composed of a 
surface-state line according to the new model (Eq.0 (dashed) 
and a quadratic background (dash-dotted) to account for tun- 
neling into bulk states. The residual of the fit is shown as a 
thin dotted line, (b) Dispersion and density of states (DOS) of 
the Lu-surface band, and (c) energy-dependent lifetime width 
resulting from the fit analysis. 



tional feature, thus this procedure is questionable. The 
model curves in Fig. [31 show the fit result without such a 
Gaussian function to emphasize the deviations. 

An inspection of the results of the band-structure cal- 
culation [Fig.^b)] explains why the above model, which 
assumes a monotonic downward dispersion of the surface 
band, could not be successful, because the Lu surface 
band follows a more complicated dispersion. The band 
maximum at fcy ^ leads to a singularity in the DOS 
that must be taken into account, and also the energy 
dependence of the tunneling probability is different for 
such an M-shaped surface band. Therefore, we introduce 
a more sophisticated model to describe the lineshape of 
the STS data. 

The simplest extension of the former model is to as- 
sume a parabolic dispersion of 4 t/l order, E(k^ ) = Eq + 
afe| + &fc|, with a > and b < 0; with this ansatz, 
E(k\\ — 0) — Eq is the local band minimum at T. The 
band maximum is located at fc max = y/ —a/26 and has the 
energy i? max = E + AE, where AE = -a 2 /4b. Thus, a 
and b can the exchanged by the more descriptive param- 
eters A; max and AE: 

E(k {[ ) =E a + 2A£(fc|,/fc max ) 2 - A£(fc||/fc max ) 4 . (4) 



For fey | ^ fcmax) this results in 
I fell (E) I = fcmax [ l : • 



En — E 



1/2 



AE 



(5) 



and the DOS of the surface band can be calculated as: 



this possibility. («) The linewidth increases with tem- 
perature, and thereby the surface state protrudes over 
the Fermi edge and becomes partially occupied. Our 
analysis of the spectrum at 60 K shows an increase of 
the linewidth by « 20 meV. In first approximation, the 
linewidth increases linearly with temperature due to en- 
hanced electron-phonon scattering [25t |27| . By extrap- 
olation, the linewidth at room temperature is estimated 
to be at least 0.2 eV, i.e. mechanism (ii) explains the dif- 
ferent results of this study and of Ref.0 as a temperature 
effect. Furthermore, taking into account the IPE energy 
resolution of w 0.2 eV, it is clear that no statement on 
the exact position of the band maximum could have been 
made in Ref. |fj. The comparison with our results points 
out the advantage of STS, which simultaneously probes 
both occupied and unoccupied states around Ep. 

Fig. |3 shows a highly resolved tunneling spectrum of 
Lu(0001) plotted together with the least-squares fit re- 
sults. In this case, the residual illustrates that the fit 
model does not sufficiently well describe the experimental 
data. As pointed out recently jig, the fit model works 
only satisfactorily for Lu(0001), if an additional Gaus- 
sian function is added to describe both the broad tails 
and the sharp peak in the spectrum. At present there is, 
however, no clear physical explanation for such an addi- 
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As the dispersion in this model is not simply quadratic, 
the approximation for the transmission coefficient in 
Eq.0is not valid anymore, and the general formula must 
be used. 

For Eq < E < E max , two solutions exist, i.e. for the 
same energy there are two contributions to the tunneling 
current with different momentum fej; < fc max and k^ > 
fcmax, respectively. This can be included by a sum of two 
transmission coefficients T a and Tt>. For E < Eq, T a is 
set to zero. 

Finally, as for the above model, the broadening of each 
state due to the finite lifetime is included by an integral 
over Lorentzian functions: 



§=PoJJ 



r(e) 



^n s (E)x 



(E- E ) 2 +(W 
[T a (E) + T b (E)} f(E - eV)dedE. 



(7) 



In summary, the fit routine possesses the free parame- 
ters: E Q , AE, fc max , I\ </> e ff, and z. At a first glance, this 
number seems to be too large for a reasonable fit. How- 
ever, due to the knowledge of the band dispersion from 
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Fig. ^b) , AE and fc max can be restricted to a region that 
is compatible with the calculated band structure. Also 
the last two parameters can be limited because the bar- 
rier height should not differ much from the work function 
of Lu 3.3 eV), and the typical tip-sample distance is 
between 5 and 15 A. 

For r, we take into account the energy dependence due 
to e-ph scattering (r e _ p h) according to the Debye model 
HE HE 0- and the quadratic energy dependence due 
to e-e scattering (r e _ e = 2/?[(7rfc B T) 2 + (E - E F ) 2 }) as 
described by Fermi-liquid theory pa. l35j|. r e _ p h is de- 
termined by the e — ph mass enhancement parameter A 
and the Debye energy Hljjj. r e _ e is described by the 
parameter 2(3 Further, an offset T is used as addi- 
tional parameter in order to include the possibility that 
r(Ep) > 0. Thus, the lifetime width is only considered to 
be energy-dependent in this model but not fc||-dependent. 
We point out that also a momentum dependence is to be 
expected, caused e.g. by interband e-e scattering with 
bulk bands, which is more probable when the excited 
state is further away from the center of the BZ. However, 
including these effects requires a much more complicated 
model [23, HiJ and is beyond the scope of this study. 

Fig.^Ja) displays the result of a fit analysis of the same 
Lu(0001) tunneling spectrum as shown in Fig. 01 but now 
with the new fit model. As can be seen from the residual, 
the model curve fits the experimental data perfectly. The 
results for the fit parameters are listed in Tab. [I] In order 
to test the reliability of the model, several STS spectra 
were fitted, and the reasonability of the parameters was 
checked. Although these parameters were treated as free 
parameters in the fit procedure, they were found to vary 
only slightly. Some values that deviate from the expected 
values of the calculation will be discussed in the following. 
However, the qualitative agreement is good enough to 
specify the dispersion of the surface band, cf. Fig. EJb). 

Comparison with the calculated band structure in 
Fig. mb) gives good compliance in AE, while the en- 
ergy position of the band - as discussed above - is again 
too high in the DFT calculation in GGA approximation. 
The difference in fc max is large. While theory puts it 
to one fourth, the band maximum in the fit is at about 
one half of the surface BZ, where the band structure is 
already dominated by bulk bands. Probably, the fc|| de- 
pendence of the fit is quantitatively incorrect, which is 
not unexpected since the model is based on planar tun- 
neling theory. The DOS of the M-shaped surface band 
[Fig- Bib)], however, shows that - due to the singularity 
at -Emax - states away from V have a dominant influ- 
ence on the lineshape of the tunneling spectrum. Thus, 
the Lu(0001) surface band is occupied at the center of 
the surface BZ, but it has much more unoccupied DOS 
slightly above Ep. Due to its unique dispersion, the sur- 
face band crosses the Fermi edge twice. The parameters 
(f) e s and z vary barely. A tip-sample distance of 11 A is 
realistic. The effective work function may be slightly too 
small but is still in a reasonable regime. 

For an estimation of the energy dependence of T, all 



relevant parameters have been set free in the fit proce- 
dure. The small offset T = T(E F ) = (2 ± 1) meV can be 
explained either by a small contribution of defect scatter- 
ing or by a poorer experimental resolution (maybe due 
to electronic noise). At the band maximum, the lifetime 
width is T(E + AE) = (8 ± 2) meV. The parameters 
for e-ph scattering are again acceptable: A = 1.4 ± 0.9 
is higher than the calculated value of Skriver et al. [37| , 
but also other experiments find lar ger electron-phonon 
mass enhancement factors up to 2.1 j38|. The Debye en- 
ergy, according to the fit result, is fkuo — (12 ± 4) meV, 
which is between the bulk value of 15.8 meV and the re- 
cently determined value for the surface-Dcbyc energy of 
10.4 meV 0. 

In contrast, the values describing e-e scattering are 
rather unprecise, particularly for Ep — E > lOmeV. 
The 2/3 parameter varies between and (unphysical) 
100 eV -1 for different tunneling spectra. For a clearer 
validation, the mean value of T is plotted versus the en- 
ergy [Fig. Uc)], the grey-shaded area marks the uncer- 
tainty region. In the relevant region ±10meV around 
Ep, where the surface-state peak occurs in STS, the vari- 
ation of T(E) is relatively small. Only for lower energies 
Ep — E > 10 meV, where the contribution of the surface 
state to the STS signal decreases rapidly, the determi- 
nation of the linewidth becomes untrustworthy, and the 
error bars get accordingly very large. In the spectroscop- 
ically crucial region, however, the error bars are smaller 
than 10 meV. To further test the model, the 2(3 factor 
was fixed to realistic values (between and 1 eV^ 1 , and 
the fit procedure was applied only within the relevant 
region of the STS spectra. The resulting parameters do 
not deviate significantly from those in Tab. ^ which em- 
phasizes that the lineshape is predominantly caused by 
the unusual DOS of the M-shaped surface band. 

We point out again that the difficulties in deriving 
precise values for the linewidth may be due to the fact 
that the model does not contain a fcy -dependence of T. 
For instance, the e-e scattering rate for an excited hole 
at the r point, i.e. deep within the local band gap, is 
predominantly determined by intraband relaxation pro- 



TABLE I: Parameters resulting from the fit analysis of the 
Lu(0001) STS data according to the new model. For com- 
parison, the "starting" values expected from calculations are 



also listed. 


parameters 


fit result 


"starting" values 


E (meV) 


-17±2 


90 


AE (meV) 

&max (A ) 


20.3 ±0.3 


34 


0.64 ±0.09 


0.27 


0efi (eV) 


2.11 ±0.05 


3.3 


z{k) 


10.9 ±0.5 


5.. .15 


To (meV) 


2± 1 





A 


1.4 ±0.9 


0.59 


huJoirneV) 


12 ±4 


15.8 


2/3 (eV- 1 ) 


0...100 


< 1 
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cesses, while for an excited hole with the same energy Eq 
at the decreasing part of the surface band, additional in- 
terband transitions into nearby bulk bands are possible. 
The interplay between intra- and interband scattering 
can be complicated, even for relatively simple systems as 
the noble metals [28l Isfij. Also, in case of the far more 
complex surface state of Lu(0001), only elaborate cal- 
culations will allow for a quantitative understanding of 
surface-state dynamics. Given the current state of the 
art, the new expanded model can be rated as a useful 
description for the STS data of Lu(0001). 

In summary, we have presented low-temperature 
scanning-tunneling spectroscopy data and density- 
functional calculations for the (0001) surfaces of La and 
Lu. Fit analyses of the spectra yield quantitative infor- 
mation about the surface-band dispersions and the elec- 
tronic lifetimes of the surface states. The surface state 



at the BZ center of La(0001) is clearly unoccupied, but 
crosses Ep due to a downward dispersion towards the 
BZ boundary. The surface band of Lu exhibits an M- 
shaped dispersion, with an occupied local minimum at T 
and an off-centered band maximum slightly above Ep. 
We demonstrated that STS results for Lu(0001) can be 
described well by an extended - but still rather simple 
- model. We conclude that STS in combination with 
ab initio calculations enables one to obtain detailed in- 
formation with high accuracy on the surface-electronic 
structure of metals near the Fermi energy. 
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